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ABSTRACT 


The  spectral  absorption  coefficient  of  heated  air  has  been  calculated  over  the 
frequency  interval  and  temperature  range  where  electronic  transitions  in  diatomic 
molecules  are  the  dominant  mechanism.  With  the  aid  of  experimentally  determined 
constants,  the  spectra  of  six  important  band  systems  have  been  theoretically 
reconstructed.  These  spectra  and  the  corresponding  line  intensities  'nave  been  pre¬ 
served  in  the  form  of  a  magnetic  tape  atlas.  Average  absorption  coefficients  and 
their  role  in  radiative  trausfer  are  discussed.  A  method  of  calculating  group  absorp¬ 
tion  coefficients  based  on  a  sampling  technique  is  described.  Digital  computer  pro¬ 
grams  based  on  this  procedure  have  been  used  to  compute  optical  transmissions 
through  heated  air.  Typical  results  of  the  calculation  are  presented. 
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I.  INTRODUCTION 


i  * 


Electronic  transitions  in  diatomic  molecules  contribute  appreciably  to  the  absorp¬ 
tion  coefficient  of  air  at  optical  frequencies  for  air  temperatures  of  a  few  thousand 
degrees^1,  2^.  Such  transitions  are  discrete  when  they  occur  between  two  bound 
states  of  a  molecule.  Continuum  transitions  occur  from  bound  to  free  states  or  between 
two  free  states  (inverse  bremsstrah  ung).  This  paper  discusses  line  absorption  arising 
in  six  prominent  molecular  band  systems  in  heated  air  and  presents  the  results  of 
some  recent  calculations  on  those  systems. 

Basic  theoretical  expressions  for  the  line  absorption  coefficient  are  presented  in 
Section  n.  Section  in  contains  a  brief  discussion  of  the  molecular  spectra  for  the 
six  band  systems:  the  Schumann- Runge  bands  of  Oz  ,  first  and  second  positive 
systems  of  ,  beta  and  gamma  bands  of  NO ,  and  the  first  negative  system  of 
N*  .  Magnetic  tape  atlases  for  the  above  spectra  are  also  described  in  this  section. 

Some  average  absorption  coefficients  and  their  utility  in  radiation  transport  will 
be  discussed  in  Section  IV,  which  also  contains  <:  description  of  an  average  transmis¬ 
sion  calculation  and  presents  sample  results  of  that  calculation. 

The  spectroscopic  notation  used  herein  is  that  of  Herzberg^,  to  which  the  reader 
is  referred  for  any  necessary  background  material. 

n.  THE  SPECTRAL  ABSORPTION  COEFFICIENT  OF  A  DIATOMIC  MOLECULE 

A  theoretical  expression  for  the  spectral  absorption  coefficient  will  be  obtained 
which  is  valid  for  electronic  transitions  in  diatom;  ?  molecules  when  the  Born- 
Oppenheimer^  approximation  may  be  applied  with  sufficient  accuracy.  Consider  an 
electronic  transition  between  two  nondegenerate  levels  or  sub-levels,  where  a  sub-level 


0 

*  « 
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is  one  of  the  components  of  a  degenerate  electronic  state.  The  absorption  coefficient 
for  such  a  transition  between  a  lower  state  L  and  an  upper  state  U  is: 

M  =  Vlr-u^Ir-u^  <!) 

where  u,  .  is  the  line  frequency,  NT  is  the  particle  density  of  the  lower  state, 

Bt  is  the  Einstein  coefficient  for  absorption,  and  F(  p )  is  a  line  shape  factor 
such  that  Jv(v)  dv  -  1  .  By  spectroscopic  convention^  the  lower  level  is  given 
by  (n",  v",  J" )  ,  where  n"  is  the  electronic  stale  indicator,  v"  is  the  vibrational 
quantum  number,  and  J"  is  the  rotational  quantum  number.  The  upper  level  is 
then  specified  by  (n1,  v',  J' )  .  Primes  always  refer  to  upper  levels  and  double 
primes  to  lower  levels.  When  v  is  expressed  in  "wave  numbers"  (cm  1)  ,  the 
Einstein  coefficient  for  absorption  is  given  by 


B 


n"  ,  v"  ,  J"  ,  d'  ,  v* ,  J' 


8tt  M"  ,  M' 


R 


M",  M'i 


3h2c 


2J"  +  1 


(2) 


where 


-*M"  ,  M' 

r\ 


R  ip* 

v",  J",M"  n' 


v»,J\M 


dr 


(») 


is  the  matrix  element  of  the  electric  moment  operator  R  =  Rg  +  Rr  of  the  electrons 
and  nuclei.  Here  M"  and  M'  are  azimuthal  quantum  numbers  numbering  the 
spatially  degenerate  rotational  levels  of  the  lower  and  upper  states,  and  the  summation 
is  over  all  possible  combinations  of  the  rotational  sub-levels  of  the  lower  with  those  of 
the  upper  state. 
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Since  the  Born-Oppenheimer  approximation  is  assumed,  the  molecular  wave 


Equation  (2)  becomes 


n  8tt  I  1 2 

Bn",n',v*,,v',J",J*  ‘  3h2c  |Re|  q(v'>v")  2J"  +  l  (8) 


The  particle  density  In  the  lower  level  is  given  by 


Nn",v",J"  ” 


N(2J"  +  l)wT exp 

n", v", J"| 

kT  J 

Q 


(9) 


where  N  is  the  total  particle  density  for  all  levels,  E  „  „  T„  is  the  energy 

n  i  v  * 

of  the  level  and  its  nuclear  spin  statistical  weight.  Q  is  the  total  partition 
function  and  may  be  written  as  the  sum  of  contributions  from  all  the  electronic 
states: 

Q  =  IX 


(10) 


Taking  into  account  electronic  (orbital  +  spin),  vibrational,  rotational,  and  nuclear 
(spin  only)  degrees  of  freedom,  we  may  write 


Q  =  Q  Q  Q 

n  electronic  vib. -rot.  muclear 


(U) 


where 


Electronic  =  u’a<2S  +  e*P 

v|eXP[ 


[^1 


*  vib- rot. 


G  (v  )hc  .  . 

o_l_  o(y) 

kT  yrol. 


(12) 


(13) 
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q: 


(v) 

rot. 


(2J  +  1)  exp 


[ 


-Fy(J)hcl 
kT“ 


dJ  ^  r^w~ 


kT_ 

hcB 


(14) 


Q 


nuclear 


(2Ia  +  1X2^  +  1) 
a 


(15) 


In  these  formulas  vQ0  Is  the  energy  of  the  lowest  vibrational  level  of  electronic  state 
n  above  that  of  the  ground  state,  Gq(v)  is  the  vibrational  term  referred  to  the  lowest 
vibrational  level,  Fy(J)  is  the  rotational  term  for  the  vth  vibrational  level  with 
its  corresponding  rotational  constant,  and  S  is  the  total  electron  spin,  is  the 
statistical  weight  for  orbital  angular  momentum  A  about  the  internuclear  axis  and 
takes  the  value  1  if  A  =  0  and  value  2  if  A  *  0.  and  1^  are  the  nuclear  spins 
of  the  two  nuclei  and  o  is  a  symmetry  number  having  the  value  2  for  homonuclear 
molecules  and  1  for  heteronuclear  systems. 

An  approximate  formula  for  Qvl^.rot  is  more  convenient  for  calculation  and 
yet  accurate  enough  for  our  problem  was  developed  by  Bethe^  and  later  corrected  by 
B rinkle>  (9) : 


1 


*vib-rot. 


1  -  exp 


-1.4388(0  \  1.4388B 


|-1.4388u'oj 


(1  ♦  rT) 


(16) 


where  y  ,  Bethe's  correction  factor  for  anharmonicity  and  non-rigidity,  has  the  value 


y 


l 

1.4388w 

0 


/2io  x 

(  o  o 

V  w 

\  0 


+ 


B 

o 


(17) 


Here  *i  and  B  are  the  vibrational  and  rotational  constants  of  the  v  *  0 
o  o 

vibrational  level,  w  x  is  the  first  anharmonlc  constant  and  a  is  the  interaction 
oo  o 

constant  for  coupling  between  rotation  and  vibration. 
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Returning  now  to  Eq  (9),  we  multiply  numerator  and  denominator  by  Qn  as 
given  by  Eqs.  (11),  (12),  and  (15)  and  obtain  after  a  little  algebra 


Nn"v"J"  ’  NtotalQ, 


V  w1(2J»  ♦  1)  exp  |-  |G0(v")  *  Fy„( J'*)| 
total  ( ?.S  +  1  )^Qnuciear^vlb.rot. 


where  use  has  been  made  of  the  relation 


E.W  ■  |*oo  *  G0<v">  *  Kc 


now  let  Pq„  =  Qn,,/Qtotaj  denote  the  fractional  population  of  the  n"  electronic 
state  and  introduce  the  notation 

8*2l0  ur,  i2 

Hn"n' ,  v'V  ,  J"J*  *  "1ST  ( 2S  +  1  )wAQnuclear  Vv  ,  J"J’ ,  v'V1^  q<v'  •  v">Sj" 


19  3 

where  Lq  =  2.6875  *  10  particles/cm  is  Loschmidt's  number. 

We  now  obtain  the  desired  form  of  Eq.  (1)  by  substitution,  using  Eqs.  (18),  (20), 
and  (21)  . 


jj(v)  -- 


NtoUlP, 

Lo^vih- 


nM  I  /”  vHJM\ 

Hr"n' ,  v’V  ,  ,I"J’  exp  i  T  /  *  F(  ^  * 


Here  the  quantities  H  and  E  are  characteristic  of  the  isolated  molecule,  I. e. . 
do  not  depend  on  the  temperature  and  density,  fince  the  bracketed  term  in  (22)  Is 


7 


dimensionless,  H  has  dimen-ions  jm"“  ,  and  F(v)  has  dimensions  cm  (discussion 
follows),  the  absorption  coefficient  u{V)  has  the  required  dimension  cm  l 

Although  Eq.  (22)  is  valid  fcr  arbitrary  torm  factor  F(v) ,  in  our  calculations  for 
for  air  we  have  assumed  a  Lorentz  line  profile,  viz. , 

F<v)  =  n  "2  7 - - - \2  W 

°  +  V  '  Vn'.v'V, 

where  <j  is  th e  line  half-width  at  half-maximum.  The  absorption  coefficient  p(  v ) 
at  v  due  to  several  contributing  (overlapping)  lines  in  just  the  sum  of  the  separate 
line  absorption  coefficients: 

--  Zm*)  (24) 

i  1 

where  1  labels  the  lines  and  the  sum  runs  over  all  lines  whose  profiles  overlap 
significantly*  at  v  . 

•Obviously  some  suitable  criterion  must  be  chosen  to  limit  the  calcination  to  some 
reasonable  number  of  lines  per  frequency  considered. 
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HI.  LINE  FREQUENCES  AND  INTENSITIES  FOR  AIR 


The  frequency,  "H"  function,  and  "E"  function  must  now  be  found  for  each  line 
which  is  to  be  included  in  the  calculation.  Theoretical  reconstruction  of  the  spectrum 
for  each  of  the  six  electronic  systems  is  accomplished  by  utilization  of  experimental 
constants  as  cited  under  the  separate  systems  below. 

(31 

Fine  structure  (composite  lines)  due  to  spin  splitting  and  A-type  doubling'  ’  has 
been  ignored  in  the  calculation  for  the  following  reasons:  (1)  the  spin  splitting  is  not 
well  known,  (2)  actual  line  widths  and  profiles  are  not  generally  known,  (3)  A-type 
doubling  increases  with  increasing  J  and  lines  are  probably  not  well  resolved  until 
very  high  J  is  reached,  and  (4)  it  is  felt  some  compensation  is  obtained  by  choosing 
the  line  half-width  large  enough  to  include  partially  resolved  lines. 

Schematic  energy-level  diagrams  are  presented  in  Figs.  1  through  6;  these  illus¬ 
trate  levels  and  transitions  for  the  band  systems  included  in  the  discussion.  The 

(3) 

conventions  used  in  the  diagrams  are  those  of  Herzberg'  .  The  resolution  due  to 
spin  splitting  and  the  A-type  doubling  in  H  states  have  both  been  greatly  exaggerated 
for  purposes  of  clarity.  Vertical  lines  (transitions)  are  drawn  in  composite  line 
groupings. 

Intensity  factors  used  in  the  calculation  arc  in  accord  with  the  rotational  sum  rule 
with  the  exception  of  the  NO  gamma  system,  whose  composite  line  intensity  factors 
cannot  be  expressed  in  simple  form.  However,  only  two  very  weak  transitions  have 
been  omitted. 

A  few  of  the  levels  for  the  N.,  second  positive  system  are  labeled  with  the  asso¬ 
ciated  nuclear  spin  statistical  weight,  which  is  written  to  the  right  of  the  level. 


3  -  3  - 

The  B  I  -  X  Z  (Schumann-Runge)  System  of  09 
_ y _ s  _ _ £ 

3  3 

With  small  dispersion,  2  -  Z  bands  consist  of  a  single  R  and  a  single  P 
branch.  Larger  dispersion  shows  each  line  to  be  resolved  into  three  components  of 
about  the  same  intensity.  There  are,  thus,  six  main  branches.  Six  very  weak  satellite 
branches  for  which  AJ  *  AK  are  also  present.  Alternate  triplets  in  the  P  and  R 
branches  are  missing,  since  the  nuclear  spin  of  the  oxygen  atom  is  zero.  Even  K 
rotation  levels  are  missing  in  the  X  state,  and  odd  K  levels  are  missing  in  the  B 
state.  In  our  calculation,  lines  will  be  treated  as  unresolved.  Figure  1  shows  two 
gioups  of  lines  which  arc  considered  to  be  composite  lines  whose  iionl-London  factors 
are  obtained  by  simply  adding  up  the  S,’s  for  the  individual  components.  These  fac¬ 
tors  turn  out  to  be  3K"  for  the  P  branch  lines  and  3(K"  +  1)  for  the  R-branch  lines. 
Rotational  terms  of  both  Z  states  are  then  given  by 

Fy(K)  -  ByK(K  +  1)  -  DvK2(K  *  l)2  (1) 

The  rotational  constants  and  Dy  for  the  ground  state  were  calculated  from 

(3'  3- 

the  molccidar  constants  of  Herzberg'  ■  For  the  B  Z  state,  the  B  ‘s  at.d  D  ’s 

*  U  V  V 

were  obtained  from  experiment^11". 

For  the  O.,  Schumann -Bunge  system,  then,  the  Hj,.  j,‘  s  are  given  by 


Vj-  V.  K- 


iSr*  I. 


"  l-  ■  /  -<  \ 
■1h‘  3  V'Wh..,/ 


‘K”.  K’  c 


(3a) 


i,.„  ...  H  qtV.v” ) 
-the  K  .  K  e 


!  ~  q  (v '  ■ v  "  > 

K  -  *  l  if  K’  K”  +  i 

K"  if  K’  K”  -  l  (2b) 


o 


where  K" 


i  .  3  , 


The  B  ^  -A  3  I()  (First  Positive)  System  ol"  N’2 
3  3 

The  structure  of  n  -  Z  bands  is  very  complicated  (Fig.  2),  particularly  if, 

as  is  the  case  here,  the  coupling  in  the  n  state  is  near  Hund's  case  "a."  Since  there 

are  three  sub-bands,  3H  -  32+  .  3il  -  3E+  ,  3n  -  3I+  ,  ana  each  sub-band  has 
o  1  i 

nine  branches  (three  for  each  triplet  component  of  the  lower  state),  there  are,  in  all, 

27  branches.  Of  these,  wc  expect  nine  branches  to  be  relatively  strong,  ten  satellite 
Dranches  to  be  somewhat  weaker,  and  eight  to  be  weak  branches;  the  latter  are  repre¬ 
sented  by  dashed  lines  in  the  figure  We  will  omit  the  latter  eight  branches  entirely. 
The  19  remaining  branches,  along  with  their  corresponding  intensity  factors,  are 
included  in  the  calculations. 

Formulas  for  the  rotational  levels  of  a  3£  state  have  been  derived  by  Schlapp^11^. 
3  + 

Now  in  the  case  of  the  A  Su  -  state  of  N2  ,  the  spin  splitting  is  very  small, 
and.  furthermore,  is  known  only  very  approximately  for  most  of  the  vibrational  levels 
Therefore,  we  ignore  the  spin  splitting  and  use  instead  the  term  formula  for  *2  states 
given  by  Eq  (1). 

The  values  oi  By  and  for  each  of  the  vibrational  levels  considered  were 

(12) 

calculated  from  molecular  constants  given  by  Naude 

3  (13) 

Rotational  term  formulas  for  the  U  state  arc  from  Butkr  and  arc  for  any 
degree  of  uncoupling: 

Fvl(J)  Bv  [j(J  ♦  1)  -  ^  -  2Z,/3Z,|  -  I\.  (j  -  J)  <-'> 

Fv2(J)  Bv  [j(j  4  U  *  tZv'T/jJ  -  Dv  (j  ♦  J)  <-*) 

Fv3(J)  Bv  [j  <J  ♦  J)  ♦  -  27. 3/  J  -  Dv(J  -  3/2  )*  (5) 

1 . 


where 


Zj  =  Yy(Yv  -  4)  +  4/3  +  4J(J  +  1)  (G) 

Z2  •=  Yy(Yv  -  l)  -  4/9  +  2J(J  -t  l)  (!) 

and  Yy  -  Av/Bv  is  a  measure  of  the  degree  of  coupling  of  the  spin  to  the  internuclear 

axis.  For  large  rotation.  F  ..  F  ,  and  F  „  go  over  into  a  case  "b"  term  series 

with  J  -  K  +  1,  K,  and  K  -  1  respectively.  Note  that  levels  F^fO),  Fy3  (0). 

and  F  3(1)  do  not  exist  because  of  the  requirement  that  J  >  Q  -  |  A+ £j  for 

Huiid's  case  ”a."  The  B  ,  D  .  and  Y  values  for  each  vibrational  level  are 
v  v  v 

M?) 

taken  from  Budo'  Since  A-type  doubling  is  small  except  for  very  large  K,  it 

has  been  neglected  in  the  above  formulas. 

II  „  for  the  first  positive  system  is  given  by 
J 


HJ" 


IhST  -iBjSiv-.V) 


2/3  if  K”  even 
4/3  if  K”  odd 


sK., 


(») 


sK" 


K”  (P-brmnch) 
K"  +  1  (R-branch) 


The  N„  C3H  -  B  II  Second  Positive )  System 
2 _ u  .  g 

3  3 

The  structure  of  II  -  II  bands  is  simple  il  both  states  bekng  either  to 


Hund's  case  "a"  or  to  llund's  case  ”b.”  If  both  II  states  belong  to  case  ''a,"  the 

3  3  3  3 

selection  rule  A  I  -  0  allows  division  into  three  sub-!>ands  H(i  -  .  11  j  '  •, 

and  .  Should  A  tape  doublmr  disregarded.  each  sub-band  has  a  strong 

h.  a  strong  P.  and  (except  for  H  -  )  a  weak  Q  branch.  If  both  H  states 

belong  to  case  "b"  or  if  both  go  over  from  rase  "a"  to  case  "b"  with  increasing 


rotation,  the  same  six  strong  hamls  occur  For  all  values  of  K  in  ease  "a  ."  anil  tor 


large  values  of  K  in  case  "b,"  the  three  P  branches  are  close  together  and  the  three 
R  branches  are  close  together,  giving  rise  to  a  characteristic  triplet  structure. 

In  the  Ng  second  positive  system,  both  II  states  belong  to  case  "a"  for  small 
K  values.  Both  go  over  to  case  "b"  for  large  values  of  K.  Figure  3  shows  the  six 
main  branches  of  the  C  il  system  and  the  corresponding  energy  level  diagram. 

A  rapid  transition  from  case  "a"  to  case  "b"  has  been  assumed.  Since  nitrogen  nuclei 
follow  Bose  statistics  (I  -  1),  the  symmetrical  levels  (s)  will  have  the  higher  sta¬ 
tistical  weight  ojj(s)  -  (21  +  1)(I  +  1)  ,  and  the  antisymmetrical  levels  (a)  will  have 
the  lower  statistical  weight  u>j(a)  -  (21+  1)1,  i.e  ,  aij(s)  =  6  and  Wj(a)  =  3.  There¬ 
fore  the  total  weight  for  the  A  doublet  is  u>j  =  9  which  is  just  1/2  of  the  "maximum" 

2 

weight  wn(2I  +1)  =18  that  one  would  expect  for  a  heteronuclear  system  of  the  same 

type  and  for  v  ch  both  nuclei  have  I  =  1^  =  1  . 

The  Q2  and  Qg  branches  have  line  intensities  that  fall  off  as  1/J"  ,  and  hence 
the  total  Intensity  residing  in  these  branches  is  small.  Therefore,  these  branches  will 
be  omitted  in  the  present  study. 

Formulas  for  the  rotational  terms  are  given  by  Eqs.  (3)  through  (7). 

If  v,  =  +  G^(v’ )  -  G^'(v")  ,  then  the  frequencies  of  the  lines  will  be 

given  by 


V  "  V'.v' ,j"(j"+i) 


v  „  .  +  F'  (J"  +  1)  -  F"(  J") 

v",v'  n'  '  n'  ' 


VJ’>  "  Mv''lv1 


,J"(J"-,)  V.v-  *  Fn<J"‘  *>  -  Fn<J"> 


for  n  1,2,3. 


I 

The  two  lines  arising  from  the  components  of  a  A-iype  doublet  overlap  in  fre-  I 

quency  and  are  identical  except  for  the  u,'j ,  one  having  the  value  6  t  lymmetrical)  and  g 

the  other  3  (antisymmetrical).  This  allows  us  to  construct  the  H  ,,  ,  „  .,,, 

n"  ,n',v"  ,v',J" ,  J' 

Hj„  j,  for  the  whole  doublet  by  simply  adding  the  H’ s  for  the  individual  compo-  | 

nents,  viz: 


8t2L 


o  1 


J",J'  3hc  ?  •  3  J 


sA 

"9h^  '  V  ,  V  1  Re  1  S.t» 


*\j'  !Re  !2  q(v’-v")  sj"(a72  *  ifk)  <"“> 


(lib) 


'nuclear 
3„  3r 


where  9/2  is  the  value  of  Q 

n 

Honl- London  factorf  for  “n-°n  systems  under  various  coupling  conditions  have 
been  given  by  Budo  (15).  For  this  study  it  will  be  sufficiently  accurate  to  use  the 
Honl-London  factors  tor  a  1 1 1  -  ^  11  type  transition: 


,R 


J"  +  1 


(12) 


(13) 


The  final  formulas  for  H  „  ,  „  then  become 

n"  ,  n' .  v".  v' ,  J",  J' 


R  8i”L  R 
,  n  on 


9hc 


vj„  !  Re  I  q(v'.  v")  (J”  ♦  1] 


(14) 


and 


for  n  1.2,3 


"I-  V  !nc  !  - 


u 


(15) 


I 

I 

I 
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The  values  of  B  ,  D  ,  and  Y  for  the  C  n  and  B  n  states  of  N0  which 
v  v  v  u  g  2 

were  used  in  this  study  are  due  to  Budo^ 

The  B  22*  -  X  22*  (First  Negative)  System  of  N* 

_ U _ 6 _ _ £ 

2 

Since  2  states  always  belong  strictly  to  Hund's  coupling  case  "by  giving  rise 

to  the  selection  rule  AK  =  ±  1  with  AK  =  0  being  forbidden,  the  structure  of  a 
2  2 

2  -  2  transition  is  simple.  There  are  two  principal  branches,  an  R  branch  and 

a  P  branch,  each  of  which  may  be  resolved  into  three  components  according  to  the 

rule  AJ  =0,  ±  1.  For  one  of  these  (AJ  =0)  AJ  *  AK . 

The  energy  level  diagram  and  transitions  which  apply  to  first  negative  bands  are 

2  +  2  + 

shown  in  Fig.  4.  Rotation  terms  of  lower  2^  and  upper  2^  states  are  both  given 
by  Eq.  (1). 

Douglas  lists  the  values  of  By  ,  Dy  ,  and  Gq(v)  used  in  the  calculation. 
Ignoring  the  spin,  the  rotational  levels  can  be  labeled  by  K  .  Grouping  lines 
arising  from  the  same  K"  level,  we  have 


H 


K" 


K' 


16*2lo 

3hc  l'K",K'l  Re' 


2 

q 


(v\v")SK„. 


2/3  if  K"  even 
1/3  if  K"  odd 


(16) 


where 


S"  K"  +  1  if  K'  K"  +  1 
K 

K"  if  K'  K"  -  1 


The  factors  to  allow  for  K  odd  or  even  arise  because  of  the  effect  of  symmetry 
2  2 

properties  on  2-2  transitions  of  homonuclear  molecules .  Even  K"  levels 
are  symmetric  (u'j(s)  C),  while  odd  K"  levels  are  antisymmetric  (Uj(s)  3). 
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The  B  ^11  -  X  2U  (Beta)  System  of  NO 
2  2 

The  structure  of  a  II  -  n  system  is  simple  if  both  states  belong  to  Hund's 

coupling  case  "a"  or  both  belong  to  case  "b"  or  if,  as  in  the  present  cane,  both  states 

are  intex  mediate  between  cases  "a"  and  "b"  and  go  over  together  from  case  "a"  to 

case  "b"  with  increasing  nuclear  rotation.  There  are  two  sub-bands  in  the  NO 
2  2  2  2 

Beta  system:  ^ ±/2  ~  ^  1/2  ^3/2  ~  ^3/2'  Each  sub-band  consists  of  a 

strong  R ,  a  strong  P ,  and  weak  Q  branch.  Since  the  intensity  of  the  Q  branches 
falls  off  like  J  1 ,  we  will  negiect  these  branches  entirely.  Figure  5  shows  the  NO 
beta  system. 

2  2 

Both  the  B  n  and  X  n  states  of  NO  are  intermediate  between  Hund's 

coupling  cases  "a"  and  "b,"  with  the  X  state  being  nearer  case  "a"  (A  =  124.2  cm-1) 

than  the  B  state  (for  which  A  ~  32  cm  1 )  .  A  general  formula  for  the  term  values  of 
2 

a  n  state  with  coupling  intermediate  between  cases  "a"  and  "b"  has  been  derived 
by  Hill  and  Van  Vleck^17^. 


Tn(v'J)  =  Te  +  Gn(v)  +  B' 


lv[M  f- 


1  +  (- 1) #l(J) 


(n  =  1) 


T>  ■  { 
v  ' 


(17) 


(J  +  1)’  (n  =  2) 


where 


v(J) 


Y  +  Y  /4  and  Y 

V  V  V 


A/B 


(18) 


2  2 

Here  the  subscript  n  -  1  stands  fora  nj/2  n  ~  2  stands  for  ^3/2 
A  subscript  n  lias  been  affixed  to  the  vibrational  term  G(v)  since  the  different  sub¬ 
states  of  the  multiplet  may  have  slightly  different  vibrational  levels.  K  rmulas  for  the 
line  frequencies  of  the  two  R  and  two  P  branches  follow. 


1G 


Ry(J)  =  »e  +  ^  --  F|{J" 

- 

FV(J") 

(19) 

R2(j)  =  ^  +  v{l]  +  F^(J"  + 

1) 

-  Fj( J" ) 

(20) 

Pt(J)  =  +  F^(J"  - 

1) 

~  F^(J") 

(21) 

P2(J>  =  ve  +  VW  +  F£(J"  - 

1) 

-  F^(J" ) 

(22) 

where 


with  By  and  Dy  having  been  calculated  with  the  .'.id  of  molecular  constants  given  by 

Gillette  and  Eyster^18^. 

2 

As  the  H  states  approach  case  "b,"  the  selection  rule  AK  -  0,  *  1  holds.  Also 

branches  with  AK  *  A.T  are  very  weak.  Disregarding  A-type  doubling  and  neglecting 

the  weak  Q  branches  and  satellite  branches,  we  see  that  the  band  tructure  is  similar 
2 

to  that  of  a  2  -  “2  transition,  i  e. ,  4  strong  branches,  if  the  satellite  lines  are 


r 


considered  as  part  of  a  composite  line,  but  Q  branches  are  still  ignored,  the  intensity 
factors  may  be  given  with  sufficient  accuracy  for  our  purpose  by 


R  2  R  i 

SJt,  =  SJI(  =  J"  +  1 


and  the  H  factors  are  then 


4ir2  L 


H 


1JM  ,J'  3hc  ^Re'  lJ",J'q(v',v")SJ" 


(26) 

(27) 


(28) 


The  A  2S  -  X  2  n  (Gamma)  System  of  NO 
2 

The  X  n  state  of  NO  belongs  neither  strictly  to  coupling  case  "a"  nor  to  case 

"b1;  but  to  a  transition  case  which  goes  from  case  "a"  to  case  "b"  with  increasing 

o 

rotation.  However,  the  A  £  state  belongs  strictly  to  case  "b,"  and  as  a  result 

of  this  combination  the  structure  of  the  NO  gamma  bands  is  quite  complicated.  The 
2 

n  state  has  been  discussed  previously  (see  NO  Beta  system).  Figure  6  shows  levels 

and  transitions  applicable  to  the  NO  gamma  system. 

2 

The  general  term  formula  for  the  X  n  state  is  given  in  the  preceding  discussion 

of  the  NO  Beta  system.  The  quantum  number  K  has  been  formally  extended  to  levels 

2 

having  small  rotation  in  the  case  of  the  A  £  state.  Rotational  terms  for  this  state 
are  given  by  Eq.  (1),  where  ttie  rotational  constants  were  chosen  to  agree  with 
experiment^19"20^. 

Intensity  factors  for  the  twelve  branches  are  g'ven  (apart  from  an  arbitrary  constant) 
(21) 

by  Earls'  \  These  formulas  were  normalized  to  obey  the  sum  rule. 


After  combining  branches  and  satellites  and  omitting  two  weak  transitions,  six 
composite  I»nes  serve  the  needs  of  our  calculation.  Intensity  factors  for  these 
lines  are  given  in  Table  1. 


Table  1 

INTENSITY  FACTORS 


Transition  and  Branch 


Pi(J)  = 


f 2J  +  l)2  +  <2J  +  1)  U(4J2  +  4J  +  1  -  2Y) 


R  /  n  -  '2J  +  l)2  +  (2J  -T  1  )U(4J2  4  4J  +  1  -  2Y) 
k2'j  *  "  a(J  -r  1) 

R^J)  +  RQ21<J)  =  ^~8j  1  {M  ‘  1  _  u'4j2  +  4J  +  1  -  2Y)J 
Ql<J)  +  QP2i(J)  =  B(jV  1)  |«J  +  7  +  u(4j2  +  4J  +  1  -  2Y)| 
Q2(j)  +  QR12(J)  “  (w  -  1  +  U(4J2  +  4J  +  1  -  2Y)| 

p2<J>  *  =  e^TT)  |w  +  7  '  *j<4j2  +  u  +  1  ‘  2Y >} 

«fe  -1/2 
U  --  I Y2  -  4Y  *  (2J  +  1)2| 

The  H.tl  lt’s  for  the  NO  Gamma  system  are  then 


■  ssr  '".I 


where  the  S!„  are  the  composite  intensity  factors  described  above 
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Digital  computer  codes  (for  the  IBM  7090)  have  been  constructed  to  write  the 
magnetic-tape  line  atlases  for  the  six  molecular  systems  to  be  included  in  the  calcula¬ 
tion.  For  each  system,  the  spectral  lines  and  their  "H"  and  "E"  functions  are  generated 
according  to  vibrational  and  rotational  quantum  numbers  and  labeled  by  an  identification 
vector.  These  lines  are  then  merged  and  sorted  according  to  line  frequency.  For  each 
spectral  line,  there  are  three  decimal  numbers  and  one  octal  number  (the  identification 
vector)  on  tape: 

v  a  HE 

a  a 

(frequency)  (identification  vector) 

where  o  is  formed  of  the  upper  and  lower  vibrational  and  rotational  quantum  numbers, 
the  branch  number ,  and  an  integer  from  1  to  6  which  specifies  the  molecular  Bystem  to 
which  the  line  belongs. 

The  systems,  their  ’Vey"  integer  tj.  and  the  total  number  of  lines  stored  on  tape 
for  each  system  are 


V  =  1 

°2 

Schumann-Run<{e 

13,836 

2 

N2 

First  Positive 

58,476 

3 

N2 

Second  Positive 

16.750 

4 

N2 

First  Negative 

21,306 

5 

NO 

Beta 

15,760 

6 

NO 

Gamma 

25,400 

Total  Lines 

151,528 

Table  2  contains  basic  information  pertaining  to  the  band  systems.  Note  that  the 

3 

4th  vibrational  level  is  the  highest  included  for  the  C  n  state  of  N9  .  Since  the  C 
state  is  perturbed  by  another  electronic  state  above  the  v'  =  4  level  ,  v'  was  not 
extended  to  higher  values.*  A  further  restriction  on  v'  and  v"  (as  well  as  on  J' 
and  J")  is  that  the  total  vib rations',  and  rotational  energy  of  any  level  was  not  permitted 
to  exceed  the  dissociation  energy  for  the  state. 


Table  2 

BAND  SYSTEM  PARAMETERS 


Molecule  and 
Transition 

Vr 

max 

v" 

max 

Number 
of  Bands 
Included 

f-number 

lRe!2 

(x  10‘36) 

Range 
(103  cm  *) 

o2  b  V  -  x  3r 

20 

19 

217 

0.048 

3.29 

21.5-50.0 

N2  0  3V  a3Iu 

10 

11 

66 

0.02 

3.34 

1.97-19.5 

N2C  VB3,1g 

4 

10 

34 

0.07 

4.6 

21.5-36.4 

n!  B  V  -  X  V 

2  u  g 

17 

18 

106 

0.0348 

2. 8o 

14.9-32.8 

no  b  2n  -  x  2  n 

6 

16 

40 

0.008 

0.7 

17.3-49.5 

NO  A  2i  -  X  2n 

7 

16 

55 

0.0025 

0.14 

31.0-50.0 

Equation  (1)  was  used  to  fairly  high  J  (~I00)  .  The  extent  to  which  this  formula 
fails  at  high  J  is  not  known,  but  it  is  felt  that  accuracy  is  only  fair  at  J  =  80  . 


(22) 

•Weak  bands  observed  by  Y.  Tanaka  and  A.  S.  Jursa  have  been  identified  as 
v'  *  5  bands. 
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In  addition  to  the  individual  system  line  atlases ,  a  tape  consisting  of  all  the  systems 


merged  according  to  frequency  has  been  written.  This  merged  tape,  containing  all 

151,528  lines,  will  serve  as  the  spectral  line  atlas  for  air.  By  masking  techniques, 

the  origin  of  each  line  may  be  found  from  its  identification  vector,  i.e. ,  its  species, 

lower  level,  etc.  To  compute  the  absorption  coefficient,  the  atlas  tapes  must  be  used 

in  conjunction  with  population  numbers  and  partition  functions  corresponding  to  the 

desired  air  temperature  and  density.  Population  factors  may  be  obtained  by  interpolating 
(23) 

in  the  tables  of  Gilmore'  ’ .  The  vibration- rotation  partition  functions  are  given,  with 
sufficient  accuracy  for  our  problem,  in  Table  3. 


Table  3 

PARTITION  FUNCTIONS 


Molecule  and  State  Vibration-Rotation  Partition  Functions 


V'  z. 

Qj(T)  = 

0.4834  T(1  +  0.0000149  T) 

1  -  exp  (-  2256/T) 

N2  *  3< 

W2(T)  = 

0.48468  T(1  +  0.0000174  T) 

1  -  exp  (-  2081. 2/T) 

n2  b  *ng 

Q3(T)  = 

0.42670  TO  ♦  0.0000144  T) 

1  -  exp  (-  2474/T) 

Q4(T)  = 

0.3617  T(1  +  0.0000102  T) 

l  -  exp  (-  3152. 5/T) 

NO  X  2II 


o  m  -  °-409H  T(1  ♦  0.0000119  T) 
w5'  '  '  1  -  exp  (-  2719/T) 


All  of  the  Franck -Condon  factors  used  in  the  atlas  construction  arc  from 
(5-7) 

R.  W.  Nlcholla'  \  The  electronic  f-valucs  were  taken  from  Treanor  and  Wurster 

(O^  Schumann -Rungc),  Bennett  and  Dalby^20*  (N*  First  Negative),  and  the  compilation 
(2) 

of  Meycrott  ct  ai. (remaining  band  systems). 


IV.  RADIATION  TRANSPORT  AND  MEAN  ABSORPTION  COEFFICIENTS 


One  manner  in  which  the  spectral  atlas  for  air  has  been  utilized  is  in  the 
numerical  determination  of  some  mean  absorption  coefficients  for  application  in 
radiation  transport  problems.  These  coefficients  and  the  method  of  their  deter¬ 
mination  will  now  be  discussed. 


Radiation  Transport 

The  radiative  transfer  of  energy  in  a  medium  can  be  calculated  if  one  knows 
the  intensity  of  radiation  1^  as  a  function  of  position  and  time,  direction  of  ray  and 
frequency.  Alcng  a  given  ray,  1^  changes  in  a  manner  which  is  determined  by 
emission  and  absorption  of  the  radiation  by  the  material  through  which  it  passes. 
Such  changes  are  calculated  from  the  equation  of  transfer  If  the  material  is  in 
local  thermodynamic  equilibrium,  the  equation  of  transfer  is 


dl 
_ r 

ds 


P  <B  -  I  ) 

V  V  V 


(1) 


where  s  denotes  the  length  measured  along  the  ray, 


„  •  ,  -ht-AT, 

p*  -  Pi  i  I  -  r  ) 

v  v 


(2) 


is  defined  in  terms  of  the  absorption  coefficient  *  (p  ,  T)  ami 


ty-n 


2ht-3  ,  hr  At  -1 

—  ( C  -  l  ) 


Cl) 


2-1 


is  Planck's  blackbody  distribution  function.  The  primed  coeffic'c  .*?  differ  from  the 
unprimed  ones  by  a  factor  ( 1  -  e  This  distinction  arises  because  of  the 

preserve  of  the  so-called  induced  emission.  Since  that  term  is  proportional  to  1^  ,  it 
is  convenient  to  subtract  it  from  the  absorption  and  denote  the  difference,  which 

appears  on  the  right-hand  side  of  Eq.  (1)  as  . 

The  integration  of  Eq.  (1)  is  formally  straightforward  and  leads  to 


w 


where 


$ 

r  !  My(  s*  )  da* 

s 

o 


(4) 


(5) 


If  1^  has  been  calculated  for  all  r?ys  going  through  a  given  point,  the  integral  over 
all  directions  and  frequencies 


* :  IK-irh^  *  kK-  <*> 

is  the  difference  between  absorbed  and  emitted  power  per  unit  volume  and  it  deter- 
m.nes  the  net  rate  of  heating  oi  the  material 

The  formal  simplicity  of  the  above  program  of  calculation  is  unfortunately  mis¬ 


leading  because  the  numerical  *ork  is  generally  prohibitive  The  difficulty  lies  in 
the  fact  tlvst  the  optical  properties  of  air  in  ihe  temperature  rar^ge  covered  by  this 


article  result  mainly  from  transitions  between  molecular  levels.  The  spectrum 
associated  with  the  major  band  systems  consists  of  an  enormous  number  of  lines 
and  the  absorption  coefficient  fluctuates  from  large  values  at  the  line  centers  to 
small  ones  between  the  lines .  Pecause  of  these  "windows"  the  radiation  at  some 
point  generally  comes  from  points  along  the  ray  which  are  an  appreciable  distance 
further  back.  The  frequency  dependence  of  the  distance  and  with  it  1^  fluctuates 
just  as  strongly  as  p^  .  One  can  define  an  absorption  coefficient  which  varies 
smoothly  with  v  by  forming  the  average 


v+Av 


1  f 

=  ^  J  “p 


dv 


(7) 


where  Av  is  chosen  to  be  large  compared  to  the  width  of  individual  lines  but  narrow 
enough  to  contain  only  a  few  strong  lines.  Similarly  one  can  define  a  smoothly  vary¬ 
ing  average  I  When  one  uses  these  averages  one  must,  however,  be  very  careful 
because  the  quantities  p^  and  fluctuate  so  strongly  that  it  is  obviously  a  poor 
approximation  to  replace  the  average  product  p^  by  the  product  P^  1,  It  is  for 
example,  not  correct  to  determine  1  by  direct  integration  of  Kq.  (l)wilh  p^  in 
place  of  p^  Instead,  it  is  in  principle  necessary  to  'icgratc  Kq.  (t)  at  all  fre¬ 
quencies  to  obtain  a  detailed  spcctiai  distribution  before  one  can  calculate  averages 
There  are  two  limiting  situations  where  this  enormous  amount  of  compulation 
car.  be  avoided  The  one  situation  arises  when  I  is  very  much  smaller  than  It 

-  y 

which  can  only  happen  for  a  transparent  medium,  i  e  .  if  Pjl.  (  i.  being  the  star  of 
the  radiating  region)  is  uniformly  small  compared  to  unity  In  that  ease  one  can 


neglect  the  absorption  and  one  obtains 


CO 

Q  =  -  4ir  /  n  'B  dv 
v  v 
o 


(8) 


By  introducing  the  so-called  Planck  mean  absorption  coefficient 

In' B  do 

-  ^  V  v 

MP  =  B(TT  <9> 


where 


B(T) 


J  By(T)dy  =  ^ 


(10) 


this  can  be  written  as 


Q  -  «yrT4  (U) 

In  the  opposite  extreme  of  an  opaguc  region  for  which  M’<  »  1  one  can  also  simplify 
K<1  (6)  Thc  lcn**h  1  >n  relation  is  the  site  of  the  region  where  the  optical 
properties  change  by  lees  than  say  5  percent  and  the  inequality  is  supposed  to  he 
valid  at  those  frequencies  where  B(,  contributes  significantly  to  B  For  an  opaque 
region  the  transfer  of  radiation  is  most  conveniently  treated  in  the  diffusion  approxi- 
mation.  This  approximation  follows  from  the  observation  that  1(  is  almost  equal 
to  which  leads  to  an  approximation  where  one  replaces  I  on  the  left-hand  side 


I 

I 

I 


of  Eq.  (1)  by  so  that 


(12) 


dB 


dB 


v  dT 


The  derivative  =  —pjr  adds  a  small  anisotropic  contribution  to  the 
isotropic  B^  so  that  the  energy  flow  in  different  directions  does  not  completely 
cancel  out.  There  is  therefore  a  flux 


F  =  -  ~  =  vfi 
3  ^R 


(13) 


where 


i 


1 


dB 

_ v 

dT 


dv 


x  dB 

/  dT 
o 


(14) 


defines  the  so-called  Rosseland  mean  absorption  coefficient  .  The  rate  of 

K 

radiative  heating  new  becomes 


q  =  -  V  •  F  =  |v  '  VT4j  (15) 


\/ith  an  absorption  coefficient  which  varies  between  very  large  and  very  small  values 
neither  the  averaging  procedure  of  Eq.  (9)  nor  that  of  Eq.  (14)  is  justified.  It  is  there¬ 
fore  necessary  to  find  other  methods  for  reducing  the  labor  to  a  reasonable  level. 
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The  Sampling  Method 

The  average  intensity  of  radiation  in  high  temperature  air  appears  to  change 
quite  slowly  with  wave  length.  This  suggests  that  it  may  be  sufficient  tc  solve  the 
transport  equation  in  a  few  sample  intervals  of  the  spectrum  whose  widths  are  large 
compared  to  the  width  of  individual  lines  but  small  enough  to  contain  only  a  few 
strong  lines .  After  performing  the  transport  calculations  to  find  ly  at  a  large 
number  n  (say  100)  of  evenly  spaced  frequencies  p.  in  that  interval  one  can  use 
these  to  obtain  a  meaningful  average  1^  for  the  sample  interval.  Calcu¬ 

lations  of  this  type  form  an  important  part  of  this  paper  and  are  discussed  in  a 
later  subsection. 


Grouping  of  Intensities  (Part  1) 

One  can  imagine  that  the  frequencies  ^  in  a  small  interval  from  v  to  v  +  Av 
have  been  arranged  into  groups  a,  b,  c,  etc.  where  the  absorption  coefficients  of  group  a 
are  all  larger  than  those  of  group  b,  which  are  in  turn  larger  than  in  group  c,  etc.  There 
will  now  be  less  of  a  spread  between  the  absorption  coefficients  within  each  group  than 
between  all  of  them. 

The  intensities  within  each  group  are  also  showing  much  less  spread.  This 
follows  because  the  intensity  at  some  point  along  a  ray  reflects  the  condiiicne  exist¬ 
ing  approximately  one  mean  free  path  (pj^  )  1  behind  eo  that  a  small  spread  in  p^ 
implies  that  the  points  of  origin  lie  closely  together. 

If  one  defines  averages  within  each  group,  i.e.  .  p^  ,  . .  . .  and  1^  .  i^  ... . 


it  is  now  quite  reasonable  tc  approximate  the  average  products  by  products  of  the  aver¬ 
ages  like 


etc .  It  is  therefore  reasonable  to  assume  that  the  average  intensities  obey  separate 
transport  equations  like 

dT  _  __ 

=  p'  (B  -  I  )  (17) 

ds  a  v  a'  '  ' 

If  there  are  k  groups,  the  heating  term  in  Eq.  (6)  can  be  replaced  by 

KiV^V^V-'  <•*> 

It  is  actually  a  fairly  good  approximation  to  use  only  two  groups. 

An  Average  Transmission  Calculation 

Consider  a  frequency  interval  A  v  which  is  small  enough  that  gross  radiation 
transport  features  do  not  change  appreciably  over  the  interval,  yet  is  large  enough 
to  contain-many  lines  (say  ~  200).  Let  Vj  denote  a  fixed  frequency  point  within  the 
interval  and  let  a  label  a  line  whose  maximum  lies  within  that  same  interval.  A 
t-ansmisBion  point  function  may  th^n  be  defined  for  an  isothermal,  homogeneous  slab 
of  thickness  x  : 

Tr(  y  =  exp  {  -  [Ip^)  +  p^(  v.)]  x}  (19) 


the  above  expression,  pc  is  an  absorption  coefficient  containing  continuum  contribu- 
t'ons  as  well  as  other  effects  not  Included  in  the  spectral  calculation  (i.e.,  contribu¬ 
tions  from  polyatomic  molecules). 

The  average  transmission  through  the  slab  for  a  frequency  interval  Ac  centered 
at  frequency  v  may  then  be  written 

Tr^iiO  *  -r-  f  exp  { -  J  2T  p '  (  v)  +  p  (  v)J  x  }  dv  CiO) 

"  Ak  «  “  c 
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Thus  far,  results  have  been  obtained  only  for  the  case  where  pc  =  0  ,  and  so  we 
now  omit  it  from  the  equations.  An  approximation  to  the  integral  is  made  by  evalu¬ 
ating  the  argument  at  many  points  in  A  v  and  utilizing  the  trapezoidal  integration 
scheme . 

An  average  absorption  coefficient  may  now  be  defined  through  the  relation 

e'M(")x  =  frAl;(  v)  (21) 

Deleting  superscripts  we  obtain 

log  T r(v) 

H(v)  - - - -  (22) 

The  sum  over  a  may  include  line  contributions  from  eeveral  molecular  species. 

Only  in  the  optically  thin  case  can  an  average  absorption  coefficient  for  air  be 

obtained  by  summing  average  absorption  coefficient  contributions  from  the  individual 

species.  The  present  calculation  is  not  restricted  to  any  particular  opacity  region. 

Based  on  the  preceding  equations,  a  digital  computer  code  was  written  to  calculate 

the  average  transmission  of  optical  radiation  through  an  isothermal,  homogeneous 

slab  of  heated  air.  Variable  parameters  in  the  calculation  are  the  air  composition 

and  electronic  population,  slab  temperature  T  ,  l'.ne  half  width  a  ,  basic  frequency 

intervals,  and  the  number  of  frequency  points  on  the  l.ne  wings. 

Results  have  been  obtained  for  temperatures  frem  )  ,0O0“K  to  12,000"K  and 

-4 

densities  from  a-mospheric  normal  to  10  that  of  normal.  Thus  far  a  value  of 
1  cm  1  for  a  h»8  been  uBed  ’onsistently  for  all  calculations  with  the  exception  of 
a  parameter  study  tor  oxygen. 
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Intervals  of  100  cm 


density  case  =  10 


-1 


lated . 


1  size  were  picked  at  2,000  cm  1  spacing  except  for  the 
where  double  the  usual  number  of  intervals  was  calcu- 


Sample  curves  from  the  results  are  presented  in  FigB.  7-01.  Smooth  curves 

have  been  drawn  through  the  computed  points,  which  fall  at  2, 000  cm  1  intervals 

from  49.500  cm  1  on  down.  (Except  the  case  where  =  10  1  ,  as  was  pre- 

Po 

viously  mentioned.) 


Grouping  of  Intensities  (Part  2) 

One  way  of  calculating  two  group  absorption  coefficients  utilizes  the  results  on 
transmission  which  are  obtained  by  means  of  the  sampling  method.  After  obtaining 
the  transmission  through  a  plane  slab, one  can  fit  itr  dependence  on  the  slab  thickness 
x  by  the  expression 


Tr(x) 


1 

2 


(23) 


To  show  the  accuracy  of  this  fit,  we  compare  the  x-depcndence  of  the  effective  absorp¬ 
tion  coefficient 


lc*eTr 


(24) 


as  calculated  from  the  slab  calculation  and  from  the  formula  of  Eq.  (23). 
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Table  4.  Comparison  of  calculated  and  fitted  values  of  the  effective 
absorption  coefficient  for  a  typical  case  T  =  5,000‘K  , 
v  =  25,500  cm-1,  p/pQ  =  1.0 


x(cm ) 

4 

8 

16 

32 

64 

128 

256 

512 

1024 

1°5'1calc 

685 

683 

678 

670 

653 

622 

568 

488 

395 

!°5  Hflt 

681 

680 

676 

670 

653 

622 

569 

495 

434 

The  differences  at  large  values  of  x  are  not  important  because  the  transmission 
is  negligible  in  that  case. 

Some  tables  of  pa  and  for  various  cases  follow  the  sample  curves  from 
the  average  transmission  calculation. 

Continuum  Contributions 

Inclusion  of  the  continuum  in  the  calculations  awaits  a  more  accurate  analysis 
of  the  various  effects  that  are  involved.  For  the  readers  convenience,  however,  a 
means  of  including  any  desired  contribution  will  be  indicated. 

Returning  to  Eq.  (20)  and  assuming  pc  to  De  a  slowly  varying  function  of  fre¬ 
quency  over  the  interval  Av  ,  we  obtain  as  an  approximation 

_av  -»cw*  r 

Tr  =  — ^ —  J  exp{-[  x)  (25) 

Av 

By  applying  Eq.  (22)  the  total  average  absorption  coefficient  is  given  by 

n(y)  ^  Mc(  v)  +  i^(  v)  (26) 
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|  where  f  labels  that  pan  of  the  absorption  coefficient  arising  from  the  lines.  The 

I  same  procedure  can  be  followed  when  one  is  using  the  two  group  absorption  coeffi¬ 

cients  .  The  corresponding  equations  are 

^  ^  +  ^af  <27> 

and 

^  -  ^  +  ^b£  (28) 

where  the  label  I  has  been  affixed  to  those  coefficients  that  were  calculated  when 
continuum  processes  were  omitted. 

Pho'odetachment  of  0  ,  inverse  bremsstrahlung.  and  the  02  Schumann-Runge 
continuum  are  among  the  contributors  to  the  continuous  absorption  coefficient .  In 
addition,  important  polyatomic  line  transitions  might  be  best  handled  by  empirical 
means  as  a  part  of  the  continuum .  One  known  important  polyatomic  contributor 
which  is  receiving  special  attention  is  NOz  .  All  these  will  eventually  be  integrated 
into  the  ca!culationa 

Discussion  of  Sample  Results 

Sample  curves  of  T  vs  X  in  Figs  7 -38  exhibit  the  smooth  behavior  that 
allows  the  two-group  approximation  to  be  a  fairly  good  one  As  the  temperature  is 
raised  and  density  lowered,  however,  some  of  the  complicated  frequency  dependence 
^  shows  up  in  the  crossing  of  these  curves,  for  which  frequency  is  a  parameter. 

Figures  39-70  present  the  same  sort  of  results  in  a  different  fashion,  i.c. ,  in 

I 

the  form  of  an  average  absorption  coefficient. 

! 

I 

I 
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Influenced  by  the  sensitivity  of  the  fractional  species  concentrations  and  elec¬ 
tronic  level  populations  to  changes  in  temperature  (and  density), the  frequency 
dependence  of  the  transmission  becomes  very  complicated  above  3,000°K  at  densi¬ 
ties  10  1  normal  atmospheric  and  below.  The  predominant  absorption  system  at 
temperatures  from  1.000  to  3,000* K  is  the  O,,  Schumann-Runge  system  which 
is  strongest  in  the  higher  frequency  regions,  i.e. .  about  30,  GOO  cm  1  on  up.  At 
3,000’  K,  NO  starts  to  contribute  appreciably  and  increases  in  importance  as  the 
temperature  increases .  A  very  important  constituent  at  the  higher  temperatures 
(6.000  to  8,000’K)  is  Ng  .  whose  first  negative  system  exhibits  a  strong  peak 
and  valley  behavior  which  is  apparent  in  the  lower  density  plots . 

More  points  were  computed  for  the  density  case  p/p^  =  10  1  .  leading  to 
lower  energy  information  not  on  the  other  graphs.  In  particular,  the  point  at 
18, 500  cm  1  shows  a  relatively  more  opaque  region.  This  is  apparently  the  high- 
energy  end  of  the  N2  first-positive  contribution  and  is  to  be  expected,  according 
to  the  work  of  Meyerott,  Sokoloff,  and  Nicholls^. 

Where  gross  energy  transfer  is  being  considered,  it  must  be  remembered 
that  the  vibration -rotation  spectrum  of  NO  ,  which  is  especially  important  at 
lower  frequencies  for  temperatures  up  to  3,000*K.  has  not  yet  been  included  in 
the  calculations,  and  therefore  is  not  contained  in  these  results. 


Table  5 

TWO-GROUP  ABSORPTION  COEFFICIENTS 
T  -  5000*  K 


Table  5  (Cont. ) 
T  =  6000°  K 


P/P0  = 


1.0 

10-1 

o 

1 

to 

10-4 

Frequency 

(cm-1) 

"a 

“b 

"a 

% 

Ma 

% 

19500 

5.07'4 

1.08'4 

9. 63-6 

2. 33~6 

3. 50-25 

2.12'26 

21500 

2.44"3 

8. 93~4 

4.27-5 

1.61'5 

7. 58_22 

7. 92_23 

2.47'8 

5.37*9 

23500 

5. 68~3 

2. 26~3 

9.48_5 

3.88~5 

i.or17 

5.08'19 

9.  23_8 

1.22'8 

25500 

1. 13~2 

4. 36-3 

1.81-4 

6.86*5 

i.io-15 

1.33-16 

2. 13'8 

1.72'9 

27500 

1.81'2 

7. 01‘3 

2.85-4 

1. 08  4 

2.66"15 

LOS'17 

1.29'8 

4. 22~9 

29500 

3. 61~2 

1. 50-2 

5.59-4 

2. 27_4 

3.44'1? 

1.64-15 

2.99'8 

1.40'8 

31500 

_2 

3.76 

1.34-2 

5.71-4 

-4 

1.94 

2.49-12 

1.49-13 

1.79'8 

6.30'9 

33500 

6. 32~2 

2.77'2 

-4 

9.63 

4. 12~4 

1.21-10 

4.97-12 

1.99*8 

8.0l'9 

355C0 

1.05"1 

4.57-2 

1.7l"3 

6. 50~4 

1.12'9 

9.92"U 

3.99'8 

1.50'8 

37500 

1.42'1 

5.08'2 

2.  n-3 

7.11-4 

1.51-8 

3.50“9 

4.32'8 

1.33'8 

39500 

2.04"1 

6.12'2 

3. 12'3 

8.62"4 

1.15'7 

1.79'8 

6.52'8 

1.58'8 

41500 

3. 15'1 

1.24'1 

4.92'3 

1. 80~3 

8. 73*7 

5.24'8 

1 . 03  7 

3. 47~8 

43500 

3.79'1 

1.47*1 

5.87~3 

2.25'3 

4. 12'6 

5. 88  7 

-7 

1.22 

5.40*8 

45500 

5.  33* 1 

2.47'1 

8.56"3 

3.78  3 

2.69"5 

3.8l'6 

-7 

1.82 

7.61'8 

47500 

5.23'1 

1.88’1 

8. 08'° 

-3 

2.79 

1.46'4 

1. 30"5 

-7 

1.69 

5. 38~8 

49500 

6.45' 1 

3.20'1 

1.04'2 

-3 

4.97 

8.41"* 

1. 10’4 

-7 

2.20 

i.or7 

3G 


I 

I 

I 

|  Table  5  (Cont.) 

T  =  7000“  K 
P/P0  = 


1.0 

io- 1 

10‘2 

10-4 

Frequency 

(cm-1) 

!  ^ 

% 

'‘b 

“a 

**a 

19500 

1.02'3 

2.74"4 

2.18-5 

7. 14'6 

2.30~6 

9.34'7 

21500 

3.55'3 

1.33-3 

6.32"5 

2.69"5 

1.02"5 

3. 19 ‘6 

2, 39~8 

5. 26_9 

23500 

7.35"3 

2.90"3 

-4 

1.49 

6.44"5 

3.12“5 

7.44'6 

7.48'® 

1.30~8 

25500 

1. 30~2 

S.06'3 

3.  ll"4 

6. 94_5 

3.51-5 

1.78*6 

2. 10~8 

1.27'9 

27500 

1. S3-2 

7. 63_3 

2.95"4 

-4 

1.26 

1 . 9 1~ 5 

8. 16~6 

1.19-8 

2. 6l”9 

29500 

3. 29~* 

1.37"2 

4.92_4 

2. lo"4 

3. 27~3 

1.57'5 

3.  24  ® 

1.34'S 

31500 

3.33~2 

1.12"2 

4. 76~4 

1.61-4 

2.57~5 

9.  25_6 

1 . 20~8 

3.549 

331C3 

5. 12~2 

2.23~2 

7.03“4 

3.00"4 

3. 29_5 

1.35'5 

4.97~9 

2.04"9 

35500 

8.30"2 

3.33'2 

1 . 22~3 

4.82"4 

6.36_5 

2.43'5 

1. 05~8 

4. 10"9 

37500 

i.oo'1 

3. 48*2 

1.37*3 

4. 53  4 

6. 33  3 

2. 03~5 

9.  68  9 

3. 00~9 

39500 

1.44”1 

3.85'2 

1.92'3 

5. 31~4 

9.23’5 

2.25~5 

1. 38  8 

3.44'9 

41500 

2.06*1 

7.97-2 

2.87*3 

1.06*3 

1. 36  4 

4.70'5 

2. oe"8 

7.  14~9 

43500 

2.41'1 

9. 13'2 

3.31*3 

1.27~3 

1.56'4 

5.71-5 

2.36'8 

8.69'9 

45500 

3.25'1 

1.45-1 

4.58‘3 

1.98-3 

2. 14  4 

9.07"5 

3.27'8 

1.36'8 

47500 

2.99'1 

1.06’1 

4. 15_3 

1.42-3 

1 .94  4 

6.46'5 

2.91 

9.87’9 

49500 

3.58'1 

1 . 73_1 

5.06'3 

2.  39  3 

2. 37~4 

i.  io-4 

3.60'8 

1.66-8 

37 


Table  5  (Cont.) 
T  =  8000*  K 


P/P0  = 


1.0 

10"1 

10-2 

10-“ 

Frequency 

(cm-1) 

% 

"a 

% 

^a 

"b 

“a  ^ 

19500 

1.95*3 

6.00~4 

4.78~5 

1.72"5 

3.42*6 

1.49”6 

21500 

4. 26~3 

1.79'3 

1.29"4 

4.85-5 

1.60”5 

4.52*6 

1.78”8  4.44”9 

23500 

9. Ol'3 

4.18'3 

3.55'4 

1.12”4 

4.985 

8.87”6 

5.56”8  8.58”9 

25500 

1.82'2 

5. 76~3 

4.45~4 

9.14”5 

3. 685 

2. Ol”6 

2.04*8  1.10”9 

27500 

2. Ol"2 

8.92'3 

3. 6l'4 

1.54*4 

2.025 

6.65'6 

l.ll”8  1.88*9 

29500 

3. Ol-2 

1.31-2 

S.42-4 

2.57'4 

3.8l”5 

1.80*5 

3. 2l”8  1.30”8 

31500 

2. 85* 2 

9.90”3 

4.12*4 

1.49~4 

1.93”5 

-6 

6.60 

33500 

3. 79~2 

.0 

1.64 

4.46~4 

1.89"4 

1.35'5 

5.68”6 

35500 

6. 47~2 

2. 61~2 

9.02’4 

3.58'4 

3. 10”5 

1.22*5 

37500 

6.86'2 

2. 29*2 

-  4 

7.92 

-4 

2.63 

2.42”5 

7.73"6 

3950C 

9.30"2 

_•> 

2.56  “ 

1.10-3 

2.88'4 

3.30'5 

8.  74  6 

41500 

1 . 3 1~ 1 

_  T 

4.92  “ 

1.53'3 

5.77”4 

4.70”5 

1.70”5 

43500 

1.47'1 

5.61'“’ 

1.74-3 

6. 52" 1 

S.26”5 

1.98'5 

45500 

1 . 9 1 _  1 

8.34'2 

2. 273 

9. 7l”4 

6.93'5 

2.91”5 

47500 

1.69"1 

6.  07~~ 

2. Ol”3 

6.95  * 

6.04”5 

2. 12'5 

49500 

1.97*1 

9.41  “ 

2. 34  3 

i.u"3 

7.13”5 

3.33"5 
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OPTICAL  PATH  LENGTH  (CM) 
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OPTICAL  PATH  LENGTH  (CM) 
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OPTICAL  PATH  LENGTH  (CM) 
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OPTICAL  PATH  LENGTH  (CM) 


OPTICAL  PATH  LENGTH  (CM) 
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fig-  17  Average  transmission  of  optical  radiation  through  a  slab  of  heated 
air  as  a  function  of  optical  pad-,  length. 
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Fig.  19  Average  transmission  of  optical  radiation  through  a  slab  of  heated 
air  as  a  function  of  optical  patn  length. 
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OPTICAL  PATH  LENGTH  (CM) 


OPTICAL  PATH  LENGTH  (CM) 
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OPTICAL  PATH  LENGTH  (CM) 


OPTICAL  PATH  LENGTH  (CM) 


Fig  30  Average  transmission  of  optical  radiation  through  a  slab  oi  heu*. 
air  as  a  function  of  optical  path  length. 
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OPTICAL  PATH  LENGTH  (CM) 
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OPTICAL  PATH  LENGTH  (CM) 
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FI*.  37  Average  transmission  of  optical  radiation  through  a  slab  of  heated 
air  as  a  function  of  optical  path  length. 
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ABSORPTION  COEFFICIENT  (CM 


FREQUENCY  ( tO*  CM*‘  ) 


Pig.  39  The  avsrage  absorption  coefficient  of  healed  air  aa  a  function  of 
frequency  for  varioua  optical  path  lengths. 
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Fig.  40  The  average  absorption  coefficient  of  heated  air  a*  a  f'^xietton  of 
frequency  for  various  optical  path  lengths. 
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Ftf  42  The  avenge  absorption  coefficient  of  healefl  air  as  a  function  of 
frequency  for  various  opt.cr.l  path  lengths. 
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Fig.  43  The  average  absorption  coefficient  oi  healed  air  as  a  fum.tion  o 
frequency  for  various  optical  path  lengths 
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Fig.  45  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig.  47  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig.  50  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig.  51  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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ABSORPTION  COEFFICIENT  (CM 


Fig.  53  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig.  54  The  average  absorption  coefficient  of  no  a  ted  air  at  a  function  of 
frequency  for  various  optical  path  lengths. 
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FI*.  5«  The  av»ra*c  aKvorption  rorffu  tens  c.f  heated  air  -»  «  function  of 
frequency  for  variou*  optical  path  lengths. 


Fig.  58  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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I  Fig.  59  The  average  absorption  coefficient  of  heated  air  as  a  function  of 

I  frequency  for  various  optical  path  lengths. 
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Fig.  60  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig.  6 1  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  foi  various  optical  path  lengths. 
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Fig.  62  The  average  absorption  coefficient  oi  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig.  84  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  Unwins. 
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F'g.  65  The  average  absolution  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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ABSORPTION  COEFFICIENT  (CM 


Fig.  66  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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ABSORPTION  COEFFICIENT  (CM- 


Fig.  67  The  average  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig.  68  The  avenge  absorption  coefficient  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths. 
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Fig-  69  The  average  ahsorplion  im;  of  heated  air  as  a  function  of 

frequency  for  various  optical  path  lengths. 
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Ft*.  70  Th«  average  absorption  ooefficirr'  of  heated  air  as  a  function  of 
frequency  for  various  optical  path  lengths 
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rl<.  71  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  as  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  is  given  in  the  figure. 


Fig.  72  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  a3  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  is  given  in  the  figure. 


Fig.  73  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  a"?  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  is  given  in  the  figure. 


Average  transmission  of  optical  radiation  through  a  slab  of 
heated  air  as  a  function  of  frequency  for  various  temperatures 
The  optical  patn  length  is  given  in  the  figure. 


Fig.  77  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  as  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  is  given  in  the  figure. 


Fig.  81  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  as  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  is  given  in  the  figure. 
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Fig.  82  Avenge  transmission  of  optical  radiation  through  a  slab  of 

heated  air  as  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  Is  given  in  the  figure. 
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Fig.  83  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  as  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  is  given  in  the  figure. 


transmission  of  optical  radiation  through  a  stab  of 
i r  as  a  function  of  frequency  for  various  temperatures. 
;al  path  length  is  given  in  the  figure. 


Fig.  87  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  as  a  function  of  frequency  for  various  temperatures 
The  optical  path  length  is  given  in  the  figure. 


Fig.  91  Average  transmission  of  optical  radiation  through  a  slab  of 

heated  air  as  a  function  of  frequency  for  various  temperatures. 
The  optical  path  length  is  given  in  the  figure. 
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